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Abstract: Cascaded Brillouin scattering is increasingly becoming of interest in many areas of pho¬ 
tonics (e.g. see [T] and [2]) and has been studied experimentally by many groups. In gain assisted 
experiments tor example up to nearly 800 distinct lines have been observed[3] while in passive cavities, 
mode locking of a few Brillouin lines has recently been demonstrated^. We show that although each 
pair of lines interact through the creation of an acoustic phonon the resulting acoustic field nevertheless 
can be written as the sum of only two waves (a forward and backward propagating wave). This new 
analysis results in a set of strongly coupled amplitude equations for the fields that are derived here 
for the first time in the context of optical fibres. More-over we show that previous authors have mis¬ 
understood the nature of the acoustic field resulting in incorrect sets of coupled equations. 

Index Terms: Brillouin Scattering, Fiber Optics 


Stimulated Brillouin scattering (SBS) in transparent materials has long been studied due to its 
ubiquity and low threshold. In optical fibres for example the threshold can be as low as a few mW 
making it the dominant nonlinear process when dealing with the propagation of narrow linewidth 
sourcesO- In Brillouin scattering an optical field generates a co-propagating acoustic field that in 
turn, changes the refractive index leading to Bragg scattering and thus a strongly reflected Stokes 
field. Conservation of energy requires that the Stokes field is downshifted by the acoustic frequency 
which in fibres is typically around 11 GHz. This reflected field can in turn be strong enough to create 
a second Stokes field leading to a cascade of multiple Stokes waves each an /c/enf/ca/frequency 
separation. Cascaded Brillouin scattering is especially strong in fibre resonators where anywhere 
up to nearly 800 distinct Brillouin lines can be seen 0 and more recently some groups have 
demonstrated that these lines can be phase-locked resulting in a 11Ghz pulse train (see (4] and 
ID). Given the importance of and growing interest in Brillouin resonators it is perhaps surprising 
that a detailed theoretical analysis of cascaded Brillouin scattering has not been performed. We 
present here a new analysis of the acoustic fields present in the fibre and show that this results in a 
new set of coupled equations for the amplitudes of the optical fields. Moreover this analysis shows 
that previous researchers have incorrectly generalised the equations for Brillouin scattering when 
trying to analyse cascaded Brillouin generation leading to flawed results. We start by summarising 
the well known physics of simulated Brillouin scattering before extending it to cascaded Brillouin 
scattering and finally discuss the differences between our work and previous analyses. 


1. Standard Analysis of Brillouin Scattering 

SBS is a three wave process in which a strong pump field P creates a co-propagating phonon field 
Q and a backwards propagating Stokes field Bq that is down shifted by the acoustic resonance 
frequency H of -- llGHz. For the case of a single Stokes wave the density fluctuations p in the 
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fibre can be written as |5] 

p{r,t) = +C.C. 


( 1 ) 


where Q is the siowiy varying acoustic fieid at frequency propagating in the +z direction with 
wave vector kA- The acoustic mode profiie is given by Fa and is assumed constant down the fibre. 
The phonon fieid Q acts to coupie the pump and Stoke’s wave ieading to the standard coupied 
ampiitude equations: 


dP 1 dP _ 
dz Vg dt 
dBo 1 dBo _ 
dz Vg dt 
dQ dQ 


—aP + ikBqQ 


—aBo + inPQ* 
-TaQ + 


( 2 ) 

(3) 

(4) 


where Vg is the group veiocity of the eiectric fieid at frequency wq, a is the ioss and k represents 
the coupiing strength. The acoustic wave propagates at a veiocity of Va and is strongiy damped by 
F/i and is driven by the product of the eiectric fieids through the eiectrostriction term represented 
by 5. Aeff represents the effective area of the fibre. 

In anaiysing Briiiouin generation it is usuai to make the assumption that the fieids are time- 
independent and furthermore the spatiai derivative of Q can be dropped due to the siow speed 
of the phonon fieid in which case Eq. H) can be soived for Q and the usuai equations for the 
power of the pump and stokes waves can be derived 0 

^ = -gslpls - ‘^alp (5) 

= gslpls - ‘^ala ( 6 ) 

dz 

where gs is the standard Briiiouin gain and Ip = PP* and = B^Bl. 

In fibres and in particuiar resonant Briiiouin iasers it is common that the Stokes wave becomes 
sufficientiy intense that it creates a second phonon wave and thus a second Stokes wave that is 
co-propagating with the pump and offset by twice the acoustic resonance frequency (2ri). Denoting 
the second acoustic wave by R and the second Stokes wave as Fq we get the coupied equations: 


dP 1 dP 
dz Vg dt 
dBo 1 dBo 
dz Vg dt 
dFp 1 dFp 
dz Vg dt 
dQ , dQ 
dt ^'^‘'dz 
dR dR 
dt dz 


—aP + inBpQ 

-aBp + iK {PQ* + FpR) 

—aFp + IK {BpR*) 

-^aQ + -^PBp 
i 5 

-TaR+-, - BpF* 

^eff 


(7) 

( 8 ) 
(9) 

( 10 ) 

( 11 ) 


Again in the time-independent case and when the spatiai dervatives of the phonon waves can be 
negiected we can soive for Q and R as 



i AAeff 

(12) 


R = ^z^^BpFp*. 
i AAeff 

(13) 
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Substituting the expressions for Q and R into the Eq. 0-(l9) and assuming continuous wave 
interactions ieads to the coupied power equations 


dip 

dz 

diso 

dz 

disi 

dz 


dsIpIsO ‘loilp 
Ob {Ipiso — Isolsi) — 2 q ;/ s 0 


93 (IsoIsi) — ‘lotisi 


(14) 

(15) 

(16) 


where Ip is the power in the pump beam, Rq is the power in the backwards propagating first 
Stokes wave and Isi is the power in the forward propagating 2nd Stokes beam. We emphasise 
here that the correct equations are Eqs. fTt-flTt rather than the coupied power equations which 
do not describe the fuii dynamics of the system (see Ref. |6] for exampie). Furthermore once a 
third stokes wave is introduced, then as we wiii show beiow, it is no ionger possibie to write down 
a set of coupied power equations and thus any anaiysis that begins with such a set of equations 
is wrong. 


2. Generation of Multiple Brillouin Lines 

Continuing the anaiysis from the previous section suppose that there are now 27V - 1 Briiiouin 
iines generated in a fibre (TV can be over 390 see for exampie Ref. O ). The fundamentai question 
is how many acoustic waves are present in the fibre? Since each pair of iines both generate and 
are coupied by an acoustic wave it might initiaiiy appear that there are 27V - 1 acoustic waves 
present of which TV are propagating in the forward direction and TV - 1 are propagating in the 
backward direction. Conservation of energy and momentum impiies that the frequency of the 
jth acoustic wave wiii be given by Stokes shift between the jth and (j + l)th eiectromagnetic 
wave. Thus we can write the totai acoustic perturbation in the fibre as 

p{z,t) = FA{x,y) I 

\ T "I 

where Qj and Rm represent the forward and backward propagating acoustic waves respectiveiy. 

In a fibre or waveguide the Briiiouin gain is reiativeiy narrow (^ 17 MHz in siiica for exampie) 
compared to the waveiength shift (11 GHz) and so the standard assumption is that aii of the 
Briiiouin iines are equally spaced In frequency impiying that 

= I = 1,2,...,2N -1 (18) 

and thus the sums in Eq. |[T7) each coiiapse to a singie term ieading to the simpiified expression 

piz,t) = FAix,y) + c.c. (19) 

for the acoustic fieid. We note that measurements of the Briiiouin shift in fibres for muitipie 
cascaded Briiiouin iines has found that Eq. |[T8) is correctQ]. 

Eq. (HH) represents the crux of our argument since if Eq. |[T8) is correct then there are only 2 
acoustic waves in the fibre: one propagating in the forward direction and one in the backwards 
direction. Looking in the iiterature we have found oniy one paper by C. Montes (7] pubiished in 1985 
that discusses this situation (in the context of piasma physics) which he caiis the strongiy coupied 
case (see Eqs. 47-50 of Ref Q). More recentiy Ogusu (8) for exampie keeps muitipie acoustic 
waves despite the fact that he assumes that each Briiiouin shift is identicai. This approach is 
foiiowed by Buttner et ai. |2] who in their suppiementary information (Eq. SI) state that the each 
Briiiouin wave is shifted by an identicai amount. A third recent exampie is Yuan et ai. (9] who 
impiicitiy adopt this modei in their numericai anaiysis of their resuits. 
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We argue here that if there is oniy a singie acoustic wave propagating in each direction in the 
fibre then it must be driven by muitipie sets of interacting fieids. This has important impiications for 
the derivation of the correct set of equations to describe cascaded Briiiouin scattering. In particuiar 
it impiies that uniike the case for 2 or 3 interacting waves coupied power equations cannot be 
used and instead coupied equations for the fieids are required and the correct equations which 
are a generaiisation of the work of Montes [7] are presented beiow. 

In line with the notation introduced above we denote the pump beam by P, the backwards 
propagating Stokes waves by Bi and the forward propagating Stokes waves by Fi where i runs 
from 0 to N. Using this notation the correct equations are: 


dP 1 dP 

dz Vg dt 

dBo 1 dBo 
dz Vg dt 


—aP + ikBqQ 
-aBo + in {PQ* + FqR) 


dB, 1 dBi 

dz Vg dt 

^ 1 ^ 

dz Vg dt 


—aBi + in {Fi_iQ* + FiR) 


—oiFi + in [B^R + Bi^iQ) 


dFN 1 dFN 
dz Vg dt 


aFjsi + in [BiR*) 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


for the optical fields (assuming an even number of fields). While the acoustic waves satisfy 


dQ dQ 

dt dz ~ 


TaQ 


TaR 


Making the standard approximations for the 



N-l \ 

PB* + 


(25) 


acoustic waves leads to the equations for Q and R: 


Q = i 


R = i 


^A^eff 

S 


N-l 

PK + E 

i-0 


N 




Em 


(27) 

(28) 


which can then be substituted back into the equations for the fields [Eq. |20j-Eq. |24)] if desired. 
In that case it can be easily shown that the coupled equations conserve the total photon flux 
through the system. 


3. Discussion and Conciusions 

The new equations Eq. |20)-l|26) are the correct and obvious generalisation of the coupled 
equations for the case of two Stokes waves when there are only two acoustic waves propagating 
in the fibre. However in the literature most authors have chosen, incorrectly, to generalised the 
coupled power equations [Eq. HD -Eq HD] or equivalently to include a separate acoustic wave 
for each new pair of Stokes waves. Not surprisingly, since this fundamental physical picture of 
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the interaction is wrong, it ieads to an incorrect set of equations and thus the resuiting anaiysis 
cannot be trusted. We emphasise here the point made by Montes in his originai articie that in 
case of more than 3 interacting opticai fieids it is not possibie to write down coupied equations 
for the evoiution of the power in the fieids and rather the coupied ampiitude equations must be 
used. 

A direct consequence of the presence of oniy two acoustic fieids is that aii opticai fieids are 
coupied together via the two phonon fieids. An effect of this is to aiter the threshoid required to see 
muitipie iines since each fieid has a growth term proportionai to the pump P (forward propagating 
fieids) or its compiex conjugate (backward propagating fieids). Thus one wouid expect that there 
wouid be more Briiiouin iines in a resonator than wouid be expected using the anaiysis of Ogusu 
for exampie. We are currentiy iooking at numericaiiy soiving the fuii coupied equations and wiii 
present a more detaiied numericai anaiysis of the difference between the two sets of equations 
in a iater paper. 

We aiso note that it is reiativeiy triviai to extend the above equations to inciude Kerr noniinear 
terms and aiso to describe ring resonators or Fabry-Perot cavities such as those treated by both 
Buttner et al. and by Ogusu. Treating a Fabry-Perot cavity requires that the number of opticai 
fieids is doubied since each Stoke wave is reflected by the cavity mirrors giving rise to a forward 
and backward propagating wave at the same frequency. However even in this case the number 
of acoustic waves remains fixed at 2 and so we do not treat this case any further. 

In conclusion we have highlighted a common error in the analysis of cascaded Brillouin genera¬ 
tion in optical fibres and presented the correct analysis based on correctly identifying the acoustic 
fields in the fibre (following earlier analysis by Montes). We hope that this work will be of use to 
other researchers working in the field. 
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